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Abstract. Wc investigate the regularizing cfTcct of adding small fractional Laplacian, with crit- 
ical fractional exponent ^, to a general first order HJB equation. Our results include some 
regularity estimates for the viscosity solutions of such perturbations, making the solutions clas- 
sically well-defined. Most importantly, we use these regularity estimates to study the vanishing 
viscosity approximation to first order HJB equations by i-Laplacian and derive an explicit rate 
convergence for the vanishing viscosity limit. 



1. Introduction 

In the pursuit of fully nonlinear partial differential equations, the Hamilton- Jacobi-Bellman type 
equations are one of the most widely, if not the most widely, studied class of equations. Among 
others, the notion of viscosity solutions has been of immense help to achieve a deeper understanding 
of fully nonlinear PDEs. It is well documented in the literature that the regularizing effect of 
adding small diffusion to first order fully nonlinear HJB equations has played a pivotal role in the 
development and understanding of viscosity solution theory. In this article we also set out to study 
a similar problem by adding a small fractional diffusion to a class of fully nonlinear first order HJB 
equations and investigate the regularizing effect and convergence properties of such approximations. 
We are interested the following initial value problem 

iut + H{t,x,u,Vu) ==0 if (t,a;) e (0,r] X R", 
|u(0,x) =ua{x) 

and its vanishing viscosity approximation 

Vm^) +e(-A)tu^ =0 if (t, x) e (0, T] x M" 
|u(0,a;) ^uo{x). 

In the above T is a positive constant, the Hamiltonian _ff is a real valued function on R x M" x R x R 
and e > is a small positive number. The precise structural assumptions on H will be detailed in 
section 2, but roughly speaking, it is a Lipschitz continuous function in all its variables and enjoys 
some monotonicity property in u. The initial data ^0(2;) is a Lipschtiz continuous function on R". 
The number s in (|1.2p is the fractional power of the diffusion operator which is supposed to be 
ranging within [1,2]. 

Among others, a rich source for equations of type (|l.ip is the area of optimal control. The 
value function of a controlled dynamical system or that of a differential game solves an equation 
of the form (jl.ll) . Also, the equations of type (|1.2p are of paramount importance due to to their 
appearance in optimal control of stochastic dynamical systems with a-stable noise. The problem 
(|1.2p is clearly a perturbation of (jl.ip . From the optimal control viewpoint, if the controlled 
deterministic dynamical system is perturbed by a small additive Levy noise then the resulting 
value function of the perturbed control problem would satisfy an equation of type (jl.2p . Our aim 
in this article is to study the stability of such perturbation and its regularizing effect on the value 
function. 
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For 5 = 2, the problem (jl.2l) becomes the classical parabolic approximation of (|l.ip and the 
classical theory for semilinear parabolic equations applies. As a result, for s=2, the Cauchy problem 
(II. 2p is well-posed and the solution u'^ is smooth (cf. [5] ) . It is also well-known that the sequence of 
functions (it'^)e>o converges locally uniformly to a function u as e J, 0, which is characterized as the 
unique viscosity solution of There are a number of methods available (cf. [BIH]) to estimate 

the rate of convergence, and one can optimally estimate the error to be of the order . 

The case s < 2 is much less classical. The operator (— A)*/^ has the following representation (cf. 

my- 

The constant C{n, s) depends only on n and s. The above integral in (II. Sp should be understood 
in the principal value sense. Clearly, in view of (jl.3p . the problem (jl.2p is non-local in nature or, 
in other words, an integro-partial differential equation. However, the notion of viscosity solution 
does make sense for such equations and the literature addressing this notion and related issues is 
fairly well developed by now. We refer to the articles [31 [HI HH [HI [12] and the references therein 
for more on this topic. The issues investigated in these papers range from standard well-posedness 
theory to more subtle questions related to regularity. 

For 1 < s < 2, the question of regularization was first answered by C. Imbert [6J. It was shown, 
under certain conditions, that the unique viscosity solution of (|1.2p is indeed of class C^'^. In other 
words, the perturbed equation (jl.2p is classically weldefined and perturbation has the same effect 
as classical parabolic regularization. In [6], the author also gives a condition on the Hamiltonian 
H under which becomes C°°. The L°°-error bound on u'^ — u is estimated to be of the order . 
The error estimate for vanishing viscosity approximation in [6] is optimal, and it was a significant 
improvement over the earlier result in [8] which was of the order . The results in [6] are new but 
mostly along the expected lines for the following reasons. The equation (jl.2p could be seen as a 
perturbation of the fractional heat equation 

ut + e(-A)4u = 0, 

and, if s > 1, this equation has similar regularizing property. Therefore, going by classical parabolic 
regularization results, it is only natural to anticipate that (II. 2p will have smooth solutions and 
results in [6] confirm this. For s < 1, the dominant derivative is of first order and featured by the 
Hamiltonian H and therefore it is not fair to expect any further regularization. In this case, it is 
well documented in the literature ( cf. 0) that equation (|1.2p will not have smooth solutions in 
general, the viscosity solutions are at best Lipschtiz continuous for Lipschtiz initial data. 

The case s = 1, as has been rightly termed, is critical. The order of the original HJB operator 
and added (nonlocal) fractional Laplace operator in (II. 2p are same, and it is a priori not clear at 
all whether there is any smoothing effect. On the intuitive level, one is more likely to think the 
opposite that there may not be any regularizing effect after all. The problem of determining the 
regularity for this critical case is much more delicate and the strategy of |6j does not apply in this 
case. It is only recently that there has been a breakthrough by L. Silvestre[T2] on this question. In 
[T2] . the author shows if H is independent of (t, x, u) then the unique viscosity solutions of (11.21) are 
indeed C^'". In other words, the viscosity solutions are regular enough to satisfy the equation in 
the classical sense. We must admit that the techniques used by L. Silvestre are fairly nonstandard 
and delicate in nature. The regularity estimate is obtained by establishing a diminish of oscillation 
lemma for the linearized version of (|1.2p . 

In this article we will concentrate on this case of critical fractional order i.e s = 1 and extend the 
results of [ni[I2]- In other words, we want to investigate the regularity of the following problem: 

ul + H{t,x,u',\/u') + e{-A)iu'' =0 ff {t, x) e {0,T] x W 
u{0, x) — uo{x). 

The contribution in this paper has two components. In the first part we extend and adapt the 
methodology of pUj to cases where the Hamiltonian H can have dependence on {t, x, u) as well as 
Vu and prove a C^'" estimate for the viscosity solution. Secondly, we use the regularity estimate to 
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estimate the error \\u'^ — u||l°°7 which comes out be of type Ce| loge|. This resuh on error estimate 
indeed estabhshes that, for critical exponent s = 1, the error estimate for vanishing viscosity 
approximation is not hnear in e. 

The rest of the paper is organized as follows: in Section 2 we state the full set of assumptions, 
relevant technical details and state the main results. We state and prove a diminish of oscillation 
lemma in Section 3. In Section 4 we use the results from Section 3 to conclude the main results. 



2. Technical framework and main results 

We begin by introducing the notations that are going to be used in the rest of this paper. By 
C, N we mean various constants depending on the data. There will be occasions where the 
constant may change from line to line but the notation is kept unchanged. The Euclidean norm on 
any R'^-type space is denoted by | • |. For any r > and x £ R", we use the notation Br{x) for the 
open ball of radius r around x. In the case when a; = 0, we simply write Br in place of Br{0) and 
define Qr — [—r, 0] x Br- For any subset Q C M x R" and for any bounded, possibly vector valued 
function w on Q, we define the following norms, 

II I / M II II \w(t, x) — w(s, y)\ 
\w\o:^ sup \w{t,x)\, \w\o,a = \w\o+ sup — — -, 

{t,x}&Q {t,x}=£{s,v} \t - sr + \x - yr 

where a € (0, 1] is a constant. The function space C^'"{Q) is the space of bounded and differentiable 
functions w such that Dt^xW = {dtw, DVxw) is Holder continuous of exponent a. This space is 
endowed with the norm 

ll^"llci-°(Q) = \w\o + \Dt,xW\o,a- 

Denote by C"^"((5) the space of all functions on Q such that |?«|o,a < oo. Also denote the set 
of all upper and lower semicontinuous functions on Q respectively by USC{Q) and LSC{Q). A 
lower index would mean polynomial growth at infinity, therefore the spaces USCp{Q) and LSCp{Q) 
contain the functions w respectively from USC{Q) and LSC{Q) satisfying the growth condition 

\w{t,x)\<c{i + \xn 

We identify the spaces USCq{Q) and LSCo{Q) respectively with USCb{Q) and LSCb{Q); is an 
index signifying boundedness. We want the initial value problem (11.11) . interpreted in the viscosity 
sense, to be well-posed and has Lipschitz continuous solutions. To this end, we list the following 
assumptions: 



(A.l) The Hamiltonian H : 
K such that 



is continuous and there is a positive constant 



||uo||lvi.=°(R") < K- 

(A.2) There exists A > such that for all {t, e R x R" x R" and u, -y G R 

H{t, X, v,p) — H{t, X, u,p) — X{v — u). 
(A. 3) There exists a constant C > such that for aU (x,p, g) e R x R" x R" and t,s eR, 

\H{t,x,u,p)-H{s,y,u,p)\ < C{l + \p\){\x-y\ + \t-s\). 
(A. 4) For every R > 0, there is a constant such that ii p,q E Bfj{0) then 

\H{t,x,u,p) - H{t,x,u,q)\ < ArIp-qI 

uniformly in {t,x,u). 

Remark. The assumptions (A.1)|[(A.4) are natural and standard, except perhaps (A.2) where the 
Hamiltonian H is assumed to be linear in 'u'. Ideally, if H is only monotonically increasing in u 
then the initial value problem (jl.ip is well-posed. However, we are interested to investigate the 
regularizing effect of ^-Laplacian on this problem and the assumption (A.2) will be necessary for 
our methodology to work. 
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We now define the notion of viscosity solution for nonlocal equations of type (11.21) . We point 
out that there could be more than one ways to formulate the definition of sub-/supersolution of the 
equation (|1.2I) . but various apparently different formulations lead to the same notion. We use the 
formulation from (H) to define the sub- and supersolutions. To this end, we introduce the following 
quantities. For k e (0, 1), let 

N 1^ /■ {u{t,x + z)-u{t,x)) 

By the representation (|1.3p . for any k, one can rewrite e{~A)i(p as 

and define viscosity solutions as follows. 

Definition 2.1 (Viscosity solution). i.) A function u £ USCb{[a,b] x M") is a viscosity sub- 
solution of p.4p if for any (p e C^'^{[a,b] x R"), whenever {t,x) e {a,b) x M" is a global 
maximum point of w — iy9 it holds that 

ipt{t,x) + H{t,x,u{t,x),\/ip{t,x)) +Il{(p) < 

for all ne (0,1). 

ii.) A function w e iS'Cb([a, 6] x M") is a viscosity supersolution of (jl.4p if for any ip G 
C^''^{[a,b] X R"), whenever {t,x) e {a,b) x M" is a global minimum point of u — (^, it 
holds that 

ipt{t,x) + H{t,x,u{t,x),\7ip(t,x)) +Il{if) +Z'"'"(m) < 
for all KG (0,1). 

in.) A function u e Cf,([a, 6] x M") is a viscosity solution of (|1.4p if it is both a sub and 
supersolution. 

The Definition 12.11 is also applicable to the case e = i.e. when the fractional diffusion term 
is absent. In this case however, the condition n S (0, 1) becomes redundant. Note that the test 
function appears in the nonlocal part of (jl.2p and this is unavoidable due to the singular nature 
of the weight function (11.31) . Some growth assumptions are needed on the sub and 

supersolutions for the nonlocal term Z'^^'^(m) to be finite; boundedness assumption is not the most 
general but sufficient for our framework. 

As usual, any classical solution is also a viscosity solution and any smooth viscosity solution is 
a classical solution. Furthermore, an equivalent definition is obtained by replacing "global max- 
imum/minimum" by "strict global maximum/ minimum" in the above definition. We may also 
assume ip = u aX the maximum/ minimum point. Next, we give an alternative (equivalent) defini- 
tion which which will be used to prove the existence of viscosity solutions via Perron's method. 

Lemma 2.1 (alternative definition). A function v e USCb{[a, &] xR") (orv € LSCbi[a, b] xR") ) is 
a viscosity subsolution (supersolution) to (|1.4p iff for every (t, x) G (a, 6) x and (j) G CQ'^([a, b] x 
R") such that v — (j) has a global maximum (minimum) at {t,x) then 

0t(i, x) + H{t, X, u{t, x), V0(i, x)) + e(-A)^ < (> 0). 

We refer to fSl for a proof. 

Remark. The Definition [3TT] is not influenced by the choice of (0, 1) as the domain of the parameter 
K. Equivalently, one can replace (0, 1) by an interval of type (0,(5) for (5 > 0. All these different 
choices for domain of n could be proven to be equivalent to alternative definition in Lemma 12.11 
However, in order for our methodology to work, we need to be able to pass to the limit k — > and 
Definition 12.11 is formulated keeping that in mind. 
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For £ e M", define i/^ : M x M" x M x M" R as 

He{t,x,r,p) = H{t,x + i,r,p) (2.1) 
and consider tfie following perturbation of (II. 4p 

+ij£(t, a;, «,£>«) +e(-A)^t; = if (t, x) G (0, T] x R". (2.2) 
We have the following continuous dependence estimate. 



Theorem 2.2 (continuous dependence). Assume (A.l)[[(Ar4) and let u, ~v e U SCh{['^,T] x M") 
respectively satisfy 

Ut + H{t,x,u,Vu) + e{-A)iu<0, (2.3) 

+ i/f(t,a;,w, Vw) + e(-A)*u > (2.4) 

in the viscosity sens^. Also assume \Du{0,x)\ + \Dv{0,x)\ < C for some constant C . Then there 
exists a constant C , depending on the data ( including T but excluding e ), such that 

v-u< (w(0, a;) - w(0, x))+ + C\e\. 

Proof. This theorem is a special case of much more general results by Jakobsen & Karlsen [5], we 
refer to this article for a detailed proof. □ 

The standard comparison principle for sub and super-solutions of (jl.2p immediately follows as 
a consequence of Theorem 12.21 if we choose ^ = 0. 



Corollary 2.3 (comparison principle). Let the assumptions (A.1)[[(A.4) be true and u,—v € 
USCb{[0,T] X R") respectively satisfy 

ut + H{t,x,u,Vu) + e{-A)^u <0, (2.5) 

vt + H{t,x,v,Vv) + e{-A)^v >0 (2.6) 

in the viscosity sense. Also assume \Du{0,x)\ + \Dv{0,x)\ < C for some constant C . In addition, 
ifu{0,x) < v{0,x) then 

u < V in Qt. 

The comparison principle ensures the uniqueness of viscosity solution for the initial value problem 
(|1.4p . The proof for existence of viscosity solution uses the standard Perron's method for viscosity 
solution framework. However, our definition of viscosity solution is slightly different compared to 
[B] and we provide detailed proof for existence. 



Theorem 2.4 (existence). Assume (A.l 



(A. 4) There exists unique viscosity solutions u € 



Cb([0,r] X R") of the initial value problem 
Proof. We begin with the claim that without loss of generality we may assume uq G C^. 

Justification: Suppose that we have proven existence of viscosity solution for initial data i.e 
Uq e Cfc2(R"). Now, if Mo e VF^'°°(R"), then there exists a sequence (ug)fe such that e C^OS.'') 
and Uq — Mo uniformly to Mo a,s k oo. Let u™{t,x) be the solution of the equation (|1.4I) with 
initial condition u"^(0,x) — m™. Then by the continuous dependence estimate 

||M"(i,x) - M''(t,x)||^o.(Q,) < C||Mo" - Mg|U^(R„). (2.7) 

Therefore the sequence (u™(t, a;)) is Cauchy in Co([0, T] x R"), and it will converge to some function 
u{t,x) G Co([0,T] X R"). We now use the stability property of viscosity solutions and conclude 
that u{t,x) is a viscosity solution of (|1.4I) with m(0, a;) — uq{x). 

Now we prove the existence of a viscosity solution for C^-initial condition. Define 

u{t,x) ^uo{x) +Ct, (2.8) 

u{t,x)=uo{x)-Ct. (2.9) 



^ This simply means that u and v are respectively the sub and supersolutions of II1.4I I 
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Invoke the assumptions (A.1)|(A.4) and choose the constant C big enough such that 



C> \H{t,x,uo{x),Vuo{x))\+e\{-^)^ua{x)\. 

Then the functions u{t, x) and u{t, x) are respectively a sub and supersolution to (|1.4|) satisfying 
the initial condition u{0,x) = uo{x). 
Define v{t, x) as 

v{t, x) — sup {w{t, x) : w < il{t, x) and w is subsolution of (jl.4p satisfying the initial condition.} 
Next, let and v* denote the upper and lower semicontinuous envelopes of v{t,x): 
v^{t,x) = linisup{u(s,2/) : (s, y) G Br{t,x) f] [0,T) x R"} 

rlO 

and v*{t,x) = — (— From the definition it is clear that 

u < w*, V* <u and v* < w*. 
The functions u, u are uniformly continuous, and hence we must have 

u = (u)* <v* <v^< {u)* = u. (2.10) 

Therefore UH.(0,a;) = w*(0,x) — uq{x). We want to show that and v* are respectively sub- 
and supersolution of (|1.4I) . This will be enough to ensure the existence, since by the comparison 
principle 

v^, < V 

and hence v^, — v* ^ v is the sought after (continuous ) viscosity solution of (|1.4p satisfying the 
initial condition v{0,x) — uq{x). 

It is relatively straightforward to show that is a subsolution and the details are as follows. 
For every {t,x) e (0,r) x M", there is a sequence {tp,Xp,Up{tp,Xp))p such that 

lim (tp,Xp,Up{tp,Xp)) = {t,x,v^{t,x)), 

where Up is a subsolution for each p e N. Now for some G C^'^ , if — has the strict global 
maximum at it, x), then there will be sequence (sp, yp)p such that Up — (j) will have global maximum 
at (sp, ?/p) and 

lim (sp,yp,Up{sp,yp)) = (t,x,v^{t,x)) . 
Furthermore, Sp > for p large enough and from the definition of subsolution we obtain 

(l^tisp, Vp) + H{sp, yp, Up{sp, yp), V(/)(sp, yp)) + e(-A)^(sp, yp) < 0. (2.11) 

Finally, we use the continuity of the equation and pass to the limit p — > oo in (j2.11l) and conclude 
that is subslolution of (|1.2I) . 

Next, we prove that v* is a supersolution. We employ the method of contradiction and assume 
that there exists {t, x) € (0, T) x M", g C^'^ satisfying v*{t, x) = (l){t, x) and w* — has a global 
minimum at (t, x) such that 

(t)t(t,x)+H{t,x,v*{t,x),V(j){t,x))+e{-A)^(j) < 0. (2.12) 

From the definition it follows that v*{t,x) < u{t,x). We claim however that v*{t,x) < u{t,x). 
Otherwise, 0(t, x) = u{t, x) = v*{t, x) and u — (f> will have a global minimum at (t, x) and 

0t(i,x) + i7(i,a;,w*(t,a;), V0(i,x)) +e(-A)50(t,a;) > 0, 

which contradicts (|2.12p 

By the continuity of and u, there are constants 71 , (5i > such that 

+ 7i<u in B5,(i,x) C (0,r) X R". 
Moreover, by (|2.12p and continuity of the equation, there exist two constants 72,(^2 > such that 
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(0 + 7)t(s, y) + H{s, y, (0 + 7)(s, y), V(0 + 7)) + A)- cj,{s , y) < (2.13) 

for all (s, y) G 3^2(1, x) and < 7 < 72. Sine v* — (j) has a strict minimum at there 
are constants 73 and < 5q < min((5i,(52) such that i>* — > 73 on dBsg{t,x). Now set 70 = 
min(7i,72,73) and define 

_ f max((/) + 70, w*) on Bs^ (t, x) n [0, T] x M" 
1 otherwise. 

Note that w is upper semicontinuous. We argue that w is a subsolution of (|1.4p . Let {s,y) G 
(0,T) X ffi." and ip G C'^'^ be test function such that ip{s,y) = w(s,y) and w — ip has strict global 
maximum at {s,y). Depending on whether w — or w — (f) + jq a.t {s,y), either — ■0 or + 70 — '0 
has a global maximum at (s, y). In the first case, the subsolution inequality for is consequence of 
I'm, being a subsolution. In the other case 

at</>(s,y) >at7A(s,y), V</)(s,y)=V^(s,2/), 1^(0) + 1«^^(0) > 1^(0) + I«^^(V^) 

and hence by (|2.13p 

tPtit, x) + H{t, X, w{t, x)ytp{t, x)) + e(-A)^ < 0. 
Therefore w is a subsolution of (|1.4p satisfying the initial condition. At the point {t,x), we have 

'w*{t, x) > max{(/)(t, x) + 70, a;)} — (j){t, a;) + 70 = x) + eg 
i.e., w(s, y) > v{s,y) for some (s, y), which contradicts the definition of w. 

□ 

It follows as a simple consequence of continuous dependence estimate that the unique viscosity 
solution of (jl.4l) is Lipschtitz continuous. 



Lemma 2.5 (Lipschitz continuity). Assume that (A.1)[[(A.4) /loW, and let u'^ e Cf, be the 



unique 

viscosity solution of (jl.4[) mi/i < e < 1. T/ien there exists a constant L, depending on the data ( 
including T but excluding e ) such that 

+ h) ~ u''{t,x)\ < L\h\ 

for all e R" and {t, x) £ [0, T) x R". 

Proof. Obviously, the function v — u'^{t, x + h) is the unique viscosity solution of (|2.2p with {£ — h 
and) initial condition w(0, x) — uq{x + h). Hence the proof follows once we apply Theorem [221 O 

Our primary aim in this article is to extend the C^'"-regularity results in |12j to the problem 
(|1.4p . To this end, it is noteworthy that the following identity holds: (cf. [T2]) 



where TZiS are the classical Reisz transforms in R". For a function u E C^'"(M"), it follows from the 
classical C" estimates for the Reisz transforms that (— A) 2 u is C". In other words, if a solution to 
(|1.4p is proven to be C^'", it will be a classical solution. A more precise mathematical formalization 
of these facts is given as the next proposition, a detailed proof of which could be found in [12] . 

Proposition 2.6. Given u £ C^'" and the integro- differential operator 

the function Cu is a C" function and C"-norm depends on ||u||(71,q and the dimension n. 

As a result, under the assumptions (A.l) (A.4)[ the nonlinear operator Xu — ut+H(t, x, u, Vii) + 
e(— A)2M also maps u to a C" function. 
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2.1. The main results. As has been mentioned a few times already, part of our main goal in 
this article is to extend the C^'"-type regularity estimate for (|1.4p . Besides, we establish that the 
solutions of (|1.4I) converges to the unique viscosity solution of (|l.ip as e — > and estimate the rate 
of convergence. 

Theorem 2.7 {C^'" — regularity). Assume 
of (ll.4p . There exists a positive constant a , depending on K, n, e and T , such that for every t > Q 
and X £ R", u is 0^'°^ at {t,x). Moreover 



(A.1)[[(A.4) and let u'^ he the unique viscosity solution 



IKII 



Ci.°((|,t]xR") 



< 



K 



\\^uo{x)\\l^). 



Theorem 2.8 (convergence rate). Assume 
respectively the unique viscosity solutions of 
depending on the data (not on e) such that 



For every e G (0,e ), let and u he 
1.41) and (|l.ip . Then there exists a constant C 



\u'-{t,x) - u{t,x)\ \ < Ce|loge| 

for all {t,x) e [0,r] X R". 

In the remaining part this section, let us outline the recipe to prove the regularity estimate in 
Theorem 12 . 71 and prove a couple of technical lemmas in connection to this recipe. Let £ G M" be a 
unit vector and assume for a moment that the solution u'^ of p.4p is smooth. Then the directional 
derivative v = diu"^ would satisfy the following linearized equation 

vt + DpH{t, x, u\Vu'').Vv + D^,H[t, x, u\ Vu)I + duH{t, x, u\SIu'-)v + e(-A)5'(; = 0. (2.14) 

In other words, v satisfies the following fractional advection-diffusion equation: 

Vt + w.Vv + \v + f{t, x) + e(-A)^t; (2.15) 

where w ~ DpH{t,x,u'^ jWu*^), f{t,x) — DxH{t, x,\/u'^) and we have used the assumption that 
H is linear in u. Keeping in line with [12,, the idea is to get a C"-type estimate for (|2.15p and 
translate that into a C^'"-type estimate of (jl.4p . 

The problem of proving Holder continuity for critical fractional advection-diffusion equations 
like (I2.15P is a delicate one, but the recent works of Cafarelli & Vasseur 1^ and Silvestre[T2] have 
contributed greatly to the understanding of this problem. In the situation when f{t,x) = A = and 
div w = 0, Cafarelli & Vasseur showed that the weak solutions to (j2.15l) becomes holder continuous 
for positive time if w G VMO. In f4], the authors use variational techniques and follow the Di- 
Giorgi type approach. In a subsequent development, Silvestre[12 uses viscosity solution approach 
and proves similar Holder continuity estimates under the only assumption that w G L°°. In this 
article we employ the later approach. 

Remember that a-priori we do not have any information on w and f(t,x) except that they are 
bounded. Also, the vector field w may not be divergence free. This makes the weak formulation by 
the means of integration by parts infeasible. However, it is possible to make sense of the inequalities 

vt-A\S/v\-B + Xv + e{-A)h<0 (2.16) 

vt+A\\7v\+B + Xv + €{~A)h>>0 (2.17) 



in the viscosity sense. If we invoke (A.l)!^^^ and select A > sup^ .j, „ [p|<|[„,||^^ \Dp{t,x,u, ■)] 



and B > C(l + ||it^||itxj) for all e, then the inequalities (j2.16l) - (l2.17l) are in perfect correspondence 
with (|2.15p . Our strategy is to follow and extend the methodology of [12] and establish Holder 
continuity for functions satisfying (|2.16p - (|2.17p and then translate it properly to C^'" estimate for 

m 

It is needless to mention that viscosity solutions are not a-priori smooth enough to undergo 
above operations. However, one can formally justify that the finite difference quotients would also 
satisfy the inequalities (|2.16p - (j2.17p . and the plan is to use this information to establish a uniform 
C"-estimate for the finite difference quotients. The next lemma makes this connection rigorous. 
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respectively satisfy 



Lemma 2.9. Assume \(ATjliAA)\ and let u,-v G USCb{[0,T] x 

ut + Hi{t, X, u, Vu) + e(-A)^u < 0, 
vt + H{t, X, V, Vv) + e{-A)iv > 
in the viscosity sense. Then {u — v) satisfies 

(u - v)t - A|V(w - w)| - B\l\ + \{u - w) + e(-A)^(u - t;) < 
in the viscosity sense. On the other hand, ifu.^v G LSCb{[0,T] x R") respectively 

Ut + Hi{t, X, u, Vm) + e{-A)iu > 0, 
Vt + H{t, X, V, Vv) + e{-A)iv < 
in the viscosity sense. Then {u — v) satisfies 

(u - v)t + A\V{u -v)\+ B\i\ + \{u -v) + e(-A)^(u - w) > 
in the viscosity sense. 



(2.18) 
(2.19) 

(2.20) 

(2.21) 
(2.22) 

(2.23) 



Before the details of the proof could be furnished, we need to introduce the notion of sup/inf 
convolution. 

Definition 2.2 (sup/inf convolutions). Given an upper semicontinuous function u{t, x) and positive 
constants (5 > 0, the sup-convolution is defined as 



u {t, x) — sup 

yeK",se[0,oo) 



1 



iy,s)~-{\x~y\' + it^.sf) 



Similarly, for a lower semicontinuous function v, the inf-convolution vs is defined as 



vs{t,x) 



inf 

j/GR'Nse[0,oc) 



1 



viy,s) + -{\x-y\^ + {t-s)^) 



(2.24) 



(2.25) 



Remark. For any upper or lower semicontinuous function on [0,T] x M", we first trivially extend 
the function on [0, oo) x K." and then define the respective sup or inf-convolution. 

We have the following lemma, the proof of which is built on ideas borrowed from [7]. 

Lemma 2.10. Let u{t, x) e USCb{[0, T] x R") and v{t, x) € LSCb{[0, T] x R") are respectively sub 
and supersolution of p.4p . Then for every z9 > 0, there exist a 5q > Q such that for all < S < Sq, 
it holds in the viscosity sense that 

(w'^)f + H{t, X, u^,Vu^) + e{-A)^u^ < i?. (2.26) 

On the other hand, for all < S < Sq , it holds in the viscosity sense that 

{vs)t + H{t, X, vs,\Ivs) + e{-A)'^vs > (2.27) 

Proof. We will provide a detailed proof for the first half i.e. (|2.26p : the proof of the second half is 
similar. 

Choose M > such that M > 2 sup |u| . For any S > 0, obviously u < u'' in [0, T] x K". Therefore 
it is easily seen that if 7 = {SoA'I)^ ^ then 

u\t, x) = sup {u{s, y) y\^ + {t~ sf) : \x - y\^ + \t - s\^ < 7} 

= sup {u{t + s, X + y) - ^ {\y\^ + {sf) : + |sp < 7} 

Clearly u{- + s,- + y) is a subsolution of 

Ut +F{t + s,x + y, u, Vu) + e(-A)3u = o (2.28) 

By (A. 2) F is linear in u and duF(t,x,u,p) is a nonnegative constant. Therefore u{- + s,- + 
— jIj/I^ + |sp is also a subsolution of (|2.28 
that u' — u{- + s,- + y) — j\y\^ + \s\^ satisfies 



y) ~ is also a subsolution of (|2.28p . Therefore, we invoke (A. 3) - (A.4) and conclude 



u't + F{t,x,u',\/u') + e{-A)^u' < C{\s\' + \y\'y^ 
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< C7 

in the viscosity sense. Choose 60, so that, C7 < i?. Then it holds in the viscosity sense it holds that 

u't + F{t, X, u', Vu') + e{-A)^u' < 1? (2.29) 

We now take the supremum and argue as in Theorem 12 .41 conclude that — sup + s,x + y) — 
+ |sp) : lyp + |sp < 7I satisfies (j2.29p in the viscosity sense. 

□ 

Proof of Lemma 12.91 By [T^, Proposition A. 5], u and vs v in the half relaxed sense as 

(5 — >■ 0. To prove the first part, it would be enough to show 

{u - v)t - A\V{u -v)\~ B\l\ + \{u - u) + e(-A)^(M -v)<2d (2.30) 

in the viscosity sense for all > 0. Hence, by the stability of viscosity solutions under half relaxed 
limit it would suffice if we prove {u^ — vs) satisfies (j2.30l) for small enough 5's. In view of Lemma 
12.101 the rest of the argument is same as [1^1 Lemma 3.2]. 

Let ip he a. test function which touches — vs at {t, x) G (0, T] x R" from above. For any S > 0, 
and —v^ are semi-convex functions, which means they have tangent paraboloid from below of 
opening j. Since the test function ip touches — vg from above at (t,x), both and —vs must 
be continuously differentiable at {t,x). Therefore, it is implied that dtv^ jdtvs^Vv^ ,\7vs are well 
defined. It is also implied that {—A)2u^ and {—A)2vs are well defined. At the point {t,x), it 
follows by direct computation that 

{u^ - vs)t - A\V{u^ ~vs)\- B\i\ + X{u^ - vs) + e(-A)^ (u^ - vs) 

<uf + Hi{t,x,u\Vu^) + e{-A)^u^ - {vs)t - H{t, x,vs,Vvs) - e{~A)hs 
<2i?, 

which clearly implies that 

ipt{t, x) - AlVifiit, x)\ - B\i\ + X{u\t, x) - vs{t, x)) + e{-A)^ip{t, x) < 2d. 
This estabUshes ((2?30l) . and thereby proves (|Z!20|) . The proof of ((2?23| is similar. 

□ 

3. The law of diminishing oscillation and C^'" estimate 

3.1. The law of diminishing oscillation. In our quest prove C^'"-type regularity, the next 
proposition plays a pivotal role. 

Proposition 3.1. Let u be an upper semicontinuous function such that u < \ in [—2,0] x R" . 
Also, u satisfies 

Ut-A|Vu|+e(-A)^u<i?o, (3.1) 

interpreted in the viscosity sense, in [—2,0] x i?2+2A- Assume further that there is a fi > such 
that 

\{u < 0}n [-2,-1] X Bi\ > fi. 

Then, for sufficiently small iSq, there is a 9 ^ (0, 1)) such that u < \ — 9 in [—1,0] x Bi. (The 
maximal value of 9 dependa on A, e and n.) 

Remark. The proof of Proposition l3.1l is a straightforward adaptation of similar results by Silvestre 
[12) . However, for the sake of completeness of our presentation, we provide the full details. 
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Proof of Proposition 13.11 Consider the following initial value ODE 

(3.2) 



^ =co|a;eBi:u(t,x)<0|-Ci7(t) for t > -2, 



i7(-2) =0, 

where 7 : [— 2, 0] i-> R be a real valued function. The solution of the initial value problem 
could be written explicitly, which is 



7(i) = J Cola; e Bi : m(s,x) < 0|e^^i(*^")ds. 



-2 

The goal is to show that, for a possibly large positive constant Ci and a possibly small positive 
constant cq, one has 

u < I ~ j{t) + 2do. (3.3) 
By Fubini's theorem, for < e [—1,0], we have 

7(0 > coe-^^^u < 0} n [-2,-1] X Bi\ > coe-^^i^. 

One can set 9 = coe~^^^^ along with the choice that < cqc^'^'^^j and the proposition follows 
from (lO). 



Choose a smooth nonincreasing function /? : K M such that I3{x) = 1 if a: < 1 and P{x) = if 
X > 2. 

Define h{t,x) = l3{\x\+At) = l3{\x\-A\t\) for t € [-2,0]. The function h{t,x) looks like a bump 
function when considered as a function of x alone. At a point x where h = (i.e. |a;| > 2 — At), 
{-A)2h < 0. Since h is smooth, (—A) is also continuous. Hence there is a constant 61 > such 
that 

(-A)5/i<0 if h<Si. 

Assume, if possible, that u{t, x) > 1 — jit) +'do{2 + t) for some point (i, x) G [—1, 0] x Also 
note that h{t,x) = 1 for aU (i,x) e [-1,0] x Bi i.e. 1 - 7(t) + i9o(2 + 1) = l--f{t)h{t,x)+'dn{2 + t) 
on [—1,0] x El. For small enough i?o, the aim is to arrive at a contradiction by looking at the 
maxima of the function 

w(t, x) = u{t, x) + i(t)h(t, x) - i9o(2 + t). (3.4) 

By our assumption, theres is a point (i,x) G [—1,0] x Bi such that w(t,x) > 1. Let (io,a;o) be the 
point where the function w achieves its maximum i.e. 

w(to,X())= max w(t,x]. 

[-2,0] xR" 

This maximum is bigger than 1, therefore it must be achieved inside the support of h. Hence 

jxol < 2 + 2|io| < 2 + 2A 
In other words, the function u satisfies p.ip at {to, xq) in the viscosity sense. Define 

ifiit, x) ^ w{ta,xa) - j{t)h{t, x) + i?o(2 + t). (3.5) 

Then {u — ip)(t, x) = w{tQ,xo) — w{t, x) and since w achieves its global maximum at {to,xo); u — ip 
has a global maximum at {tQ,xo) and u{tQ,xo) = (p{ta,xo). Therefore, for k > small enough, we 
must have 

ft{to,xo)-A\W^{to,xo)\ +r,{v){to,xo) +r'''{u){to,xo) < ^q. (3.6) 

The next task is to estimate each term on the lefthand side. To this end, notice that \WLp{to, xo)\ = 
j{t)\\7h{to,xo)\, and 

Vt{to, Xo) = -j'{to)h{to, Xq) - ^{to)ht{to,Xo) + z?o 
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- -j'{to)Hto,XQ) ~ 7(io)/3'(kl - A\t\)A + 

^ ~j'{tQ)h{to, xo) + 7(io)|V/i(io, xo)\A + do (3.7) 

The tricky part however is to obtain a refined estimate for the term 1^(95). Let us choose < k << 
1. At t = to, a function of x alone u{to^ •) + 7(to)/i(io; ■) achieves its maximum at xq. To this end, 
we denote ft = {u{to,-) < 0} n Furthermore, 



^ (i:i^)it„,x„) + r^'^iu^to, xo) 



ec{n, 1) 



B(0,k) FI"" JB(0,k)C 



(^(io,a;o + z) - </?(io,a;o)^^ _l_ f u{to, xp + z) - ^(to, a;o) 





|z|»+i 


</'(io,y) 


- </'(io,a;o) 








- 95(^0, a;o) 


\xo - 






- cp(to,a;o) 



+ / u{toy)-u{toXo)^^ 



/ ^(^o,?/)-^(to,xo) ^^ 



u{to,y) -u{to,xo) 

i TTi "2/ 



u{to,y) - u{to,xo) 
i rm "2/ 

s(a;o,«)^nn Fo - ^1"+' 

fjioiy) - '^(^O'^o) ^g^ I /" u{to,y) +j{to)h{to,y) - u{to,xo) - j{to)h{to,xo) 

\xo-y\''^'^ JBixo,K)Cnn |a;o-y|"+^ 

/, N f h{to,y) - h{to,xo) , , f u{to,y)+^{ta)h{to,y)-u{to,xo)~i{to)h{to,xo) , 

(3.8) 

Let z e ]R" such that u(toj ^) < 0, then 

m(<o, z) + j{to)h{to, z) - u{to, Xo) - -f{to)h{to, xq) < j{ta) - 1, (3.9) 

as u{to,xo) + j{to)h{to,xo) = w{to,XQ) + ?9o(2 + ^o) > 1- 
Choose Co small enough such that 7(io) < ^, then by p.9p 

u{to, z) + ■y{to)h{to, z) ~ u{to, xo) - 7(^0)^(^0, xo) < (3.10) 



2 

Therefore, for such a choice of cq, we must have by p.8p and p.lOp 



ec(n, 1) 



^ (i:(¥>)(to, 2^0) +2:^'"(")(io, 2:0) 



<^^{-A)ihito,xo) - Co\n - S(xo, ^)| (3.11) 
c(n, 1) 

where Cq is a positive constant depending on A. 

Note that as k — 0, the measure of the set \il\B{xo, k)\ — ?> We now have to consider two 
different scenarios depending on h(to,xo) and arrive at contradictions in both cases. We definitely 
have either h{tQ, xo) < Si or h{to, xo) > Si. 

In the case where h{to,xo) < Si, one has (— A)^/i(to, a^o) < 0. Hence 



(i:,{ip){to,xo)+I'-^{u){to,xo)) >ec{n,l)Co\n\B{xo,n)\ 



(3.12) 
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-l'ito)h{to,xo) + j{to)\Vh{to, xo)\A - j{to)\Vh{ta, xo)\A+I'^{ip){to,xo) +T^^"(w)(to, a^o) + ^^o < ^?o. 
In other words 

j'{ta)h{to,xo) > ec{n,l)CQ\n\B{xo,K)\, 

which is a contradiction to p.2p for any Ci if k is small enough and cq is chosen small enough so 
that it satisfies cq < ec{n, l)Co,. 

On now turn our attention to the other case h{tQ,xo) > Si. Since h{t(),xo) is a smooth function 
with compact support, we must have |(— A)2/i| < C for some C > 0. Therefore we have the 
estimate 

I^M(io,xo)+r'"(u)(to,a;o) > -eC-f{to) + co\n\B{xo, ^)\. (3.13) 
We plug (jXTO)) and dXT]) into and obtain 

-C'j{to) + co\n\B{xo, k)\ - j'{to)h{to,xo) < 0. 
We replace 7'(io) by using (|3.2[) in above and pass to limit k — > to obtain 

{Cih{to, xo) - C)j{to) + co(l - h{to, xo))\n\ < 0, 
which is contradiction under for large enough C'l as h{to,XQ) > 6i. 

□ 

Theorem 3.2 (diminishing oscillation). Let ^, C be two bounded continuous functions satisfying the 
inequalities 

6-^|Ve|+e(-A)^e<0 

Ct+A|VC|+e(-A)^C>0 

in the viscosity sense in Qi = [—1,0] x Bi. Furthermore, 

{{t, a;) e Qi : C < 0} U {{t, a;) e Qi : C > 0} = Qi. (3.14) 

Then there are universal constants 9 £ (0, 1) and ao > (depending only on A, e and n) such that 
if 

max(|C|,|CI) < 1 in Qi = [-1, 0] x Bi, 

max(|e|. Id) < 2|(4 + 4A)x|" - 1 in [-1,0] xBf 

for some < a < ao, then 

min(oscQi/(4+4A)^, oscq^ /(4+4^)C) < 2(1 - 6). 

Proof. Let i? = 4 + 4A. Consider the following rescaled versions of ^ and (: 

and <.C(-i.i)^ 

By the condition (|3.14l) 

\Bi\ 



either |{^ < 0} n ([-2, -1] x Bi)\ > 
or |{C>0}n([-2,-l]xBi)|> 



2 

\Bi\ 



2 

Without loss of generality we may assume that the former is true and it would be enough if prove 
that 
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oscQ,/^e< 2(1-0). 

The proof would be immediate if we were able to apply Proposition 13.11 to i^. Only condition 
that is missing here is that ^ needed to be bounded above by 1. To this end we define 

u ~ min(l, ^) 

and identify the inequality satisfied by u in the viscosity sense. 

Note that inside Qfli ^ ^ 1 i-e. u = ^ in Q/j. Let if be test function such that u — (p has a 
global maximum at {t,x) and u{t,x) = fit^x) where (i,a;) G Therefore ^ — (/J has a global 

maximum at (t^x) G Q2^2A- Therefore at for all k G (0, 1), 

dtip-A\V^\ + < 0. 

We are interested in estimating the quantity T^''^u(t, x) and we proceed as follows. For all k € (0, 1), 

Z^-'''l{t,x)~T-'^u{t,x) 
<^c{n,l) f (^(t,x + y)-l)^-^ 

<-Ki) / 2(|.+,r_i)^ 

<i?o, {i^o is from Proposition 13. II) 

which holds, as a result of Fatou's lemma, for small enough ao. Therefore, the function u satisfies 
the inequality 

ut - A\Vu\ + e{-A)iu < -do 

in the viscosity sense, for a suitably chosen ao, in [—2, 0] x i?2+2A. Therefore we can apply Propo- 
sition [O] to u and conclude the theorem. □ 

3.2. C^'"-regularity: the end game. We begin this subsection with the following lemma. 

Lemma 3.3. Let u be a bounded continuous function in [0,r] x M" such that it satisfies the 
inequalities 

Ut - A\Vu\ - B + \u + e(-A)5u < 0, 
ut + A\Vu\ + B + Xu + e(-A)^u > 
in the viscosity sense in [0,T] x K". Define 



C{t,x) 



[e^*u(t,a;) - f (e^* - 1) if X ^ 
[u{t,x)-Bt if A = 0. 



e^*u(i,x) + f (e^*-l) if X^O 
u{t, x)+Bt if X = 0. 



Then ^ and C are two bounded continuous functions on [0, T] x M" and ssatisfy the inequalities 

^t-A\VC\+ei~A)H<0 (3.15) 
Ct + ^|VC|+e(-A)k>0 (3.16) 

in the viscosity sense in [0,T] x K.". 
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Proof. Let {to, xq) G (0, T) x R" be a point and (p he a, test function such that (f{to, xq) — ^{to, xq) 
and S, — If has a global maximum at (to, xo). Now define 

fe-^V(<,a;) + f (f - e"^*) if A =^ 
\(p{t,x)+Bt if A = 0. 



ip{t,x) 



Then u{t,x) — 'il;{t,x) = e ^'^(^^{t,x) — ip{t,x)^, which means (^0,2^0) is also a point of global 
maximum for u — ip. Hence, 

Mh, Xo) ~ A\Wil;{to, xo)\-B + \i){to, xo) + e(-A)^(to, a;o) < 



i.e. e-^*« 



< 



- A\\/ip{to,xo)\ + e(-A) 2 ip{to, Xo) 

i.e. - A\'^ip{to,xo)\+ e{-l\)^ip{to,xo) < 0, 

which proves that (I3.15P is satisfied in the viscosity sense. The proof that (|3.16p holds in the 
viscosity sense is similar. 

□ 

Theorem 3.4. (C^ "'-regularity for the linearized problem) Let u be a bounded continuous function 
in [0, T] X R" such that it satisfies the inequalities 

ut - A\S/u\ - B + \u + e(-A)^u < 0, 

ut + A\Vu\ + B + \u + e(-A)^w > 

in the viscosity sense in (0,T) x R". Then there is a positive constant a > (depending on A,B,n 
and e) such that, for every t > and x G R" , the function u is Holder continuous of exponent a at 
{t,x). Moreover, it follows that there exist constants C and K (depending on A,n,e ) such that 

\u{t, X) - u{s, y)\ < CiMO, ■)\\l^ +K)[ ^''~ ~ (3-17) 

for allx,y&W^ and < s < t < T. 

Proof. Wc fix (to, a^o) G (0, T) and consider the following transformations of u{t, x): 

^-^^ .[u{to{t + l),xo + tox)-§{l-e-^'oit+i)-^^ if A^O 

It i .[u{toit + l),Xo + tox) - Bto{t + l)] if A = 0. 

(3.18) 

and 



at,x) 



^[u(to(i + l),xo + io.T) + f (l-e-^*«(*+i))] if A^O 

^ Au{toit^'\),xo + tox)+Bto(t+\)\ if A = 0. 

(3.19) 

where C\(^A,B ,n,t) is a positive constant, depending on the quantities in the parenthesis, to be 
chosen later and Co(B,T, X) has the following form 



CoiB,T,\) = 



- J A 



1-e-^^) if A 7^0 
BT if A = 0. 



A C"-type estimate for ^ or C at the point (0, 0) would result in a C" estimate for u at {to, xo). To 
this end, we define 

-, ei^oU+i) ^u{to{t+l),xo+tox) if A^O 

rj{t,x) = <; (ll"ll-(Q.,+2Co(B,T.A)+C.(A,B,„,.)) 

' J i Mto{t + l),XQ+tox) if A==0. 
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Then by Lemma 13.31 the functions ^ and ( are two bounded continuous functions satisfying the 
inequahties 

6-^|Ve|+e(-A)5^<0 (3.21) 
Ct + ^|VC|+e(-A)k>0 (3.22) 

in the viscosity sense in [—1,0] x R". Also note that 

max{|^|,|CI,hl}< 1- (3.23) 

Set r — j^^^- We wish to show that there exist a > such that 

msLyi(^oscQ^i^^{t,x),oscQ^^({t,x),oscQ^^ri{t,x)j < 2r"'' for ah fc^0,l,2, (3.24) 

The estimate ((3?24l) is a necessary and sufficient for Holder continuity (with exponent a) of ^, C and 
T] at (0,0). This Holder continuity of 77 at (0,0) is equivalent to Holder continuity of u at {to,xo)- 
To prove p.24[) . we follow [12^ and find a sequence (ofc, bk), fc = 0, 1, 2, such that 

flfc < ^{t,x),r]{t,x),C{t,x) < bk for all {t,x) e Qr^, 

with bk — dk — 2r"'"' where {ofe} is nondecreasing and {6fc} is nonincreasing. We employ the method 
of induction. We invoke (j3.23l) and choose ao < min(infQj ^(t, a;), infg^ C(i, x), infQ^ ri{t,x)) and 
bo > max ( supg^^ ^(t, a;), supg^ C,{t, x), supg^ ri{t, x)) such that bg — qq = 2. We now assume that 
the sequence (om, 6m) has been constructed up to some index k. We want to show the existence of 
(aA;+i,6fc+i). Define 

a(^,x) = (e(r'=x,r'=^)-^4^)r-"^ 
Ck{t,x) = {C{r''x,rH)- 

r]k{t,x) ^ {rjir''x,rh) ^ 

For any (t,x) S Qi, either ^k{t,x) < or (k{t,x) > 0. The functions ^k{t,x) and (k{t,x) 
respectively satisfy p. 211) and (I3.22p in the viscosity sense in [—1,0] x M". Also, by induction 
hypothesis, it follows that 

max{|^fc(t,x)|, |Cfc(i, a;)|} < 1 in Qi. 
Furthermore, simple calculation shows that 

max{|efc(t,x)|,|a(i,a;)|}<2r-"'"-l in Q,-^ 

for m = 1, 2, ,k. Hence 

ma^{\^k{t,x)l\Ck{t,x)\} <2\r-^x\'' -I if |.t| > 1. (3.25) 

Therefore, by Theorem 13.21 there are universal constants ao G (0, 1) and 9 e (0, 1) (depending on 
n, A, e) such that if (j3.25l) holds for any a < ao, then 

either oscg^^t < 2(1 — (?) or oscg^Cfc < 2(1 — 6*). 

In either case we have, via a simple and straightforward computation, 

/ \ 2rBT 

max[oscQ^^k{t,x),oscQXk{t,x),oscQ^r]k{t,x)] < 2(1 -9) + —— ^ .^ „ VT' 

V / [\\u\\l'^ +2C[B,T,X) + Ci(A,B,n,e)) 

where we have used the fact that r < 1. Now choose Ci{A, B, n, e) big enough such that 



2 


J 


ftfe + 


bk-^ 


2 




flfc 4 


-bk 



2rBT 

[ullioo +2C(S,r, A) + Ci(A,B,n,e)) 



Then 
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max(oscQ^£,k{t,x),oscQXk{t,x),oscQ^T]k{t,x)'^ < 2(1 - -). 
Now choose a, possibly strictly smaller than ao, such that 

Therefore we finally have 

m8i-K(^oscQXkit,x),oscQXk{t,x),oscQ^r]k{t,x)^ < 2r". 

In other words, 

max (oscQ^^^^^{t, x),oscQ^^^,C{t,x),oscQ^^^^r]{t, x)j < 2r"('=+i). (3.26) 
Therefore the pair {ak+i,bk+i) could also be chosen. □ 



With the Theorem 13.41 and Lemma 12.91 at our disposal, we can now follow the line of argument 
by SilvestrefT^ and prove Theorem 12.71 

Proof of Theorem YZ .7\ Let t e M" be a unit vector and /i be a nonzero constant. The difference 
quotient of along t at (i, x) is defined as 

u^{t,x + M) - u{t,x) 
dh^iu\t,x) = . 

For fixed h and the function dh,iu'^{x) is continuous and bounded by ||Vm'^||loo, which is bounded 
above by a constant independent of e. We now recall Lemma 1^791 and see that u = dh.iu'^{x) satisfies 
all the hypothesis of Theorem 13.41 Therefore, there exists a positive constant a such that 



WOhju (i,a:)||co.c([|^t]xR„) < — 

uniformly in h and (.. We now apply Arzela-Ascoli's theorem and and pass to the limit /i — >■ and 
conclude that diu'^{x) exists and 

C{\\\/uo\\l^+K) 



l|5iw'(t,x)||c.o,»([|,f]xR") 



< 



for all unit vectors £ £ R". In other words 

liv ^f, M\ ^ C{\\\/uo\\l^+K) 
||Vw (i,a;)||c...([|,t]xK") < -j^ • 

With this information at hand , we see that u'^{t,x) is C^'" in space for any fixed t > 0. Therefore 
H{t, X, u, Vu) + e(— A) 2 u'^ is bounded. Hence dtu'^ is bounded for t > i.e. the difference quotients 
in t 

u{t + h,x) - u{t,x) 

dhu{t,x) = 

h 

is bounded independently of h. Also, it is easy see that dhu{t,x) satisfies all the hypothesis of 
Theorem 13.41 . Therefore dhu{t,x) is bounded in C"'" uniformly in h and consequently dtu(t,x) is 
C°'" with the following estimate 

\\OtU ||c°([|.t]xR") S — ■ 

□ 
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4. Error estimate for vanishing viscosity approximation 



Proof of Theorem \2M The proof uses the doubHng of variables technique, wellknown in the viscos- 
ity solution theory. 
Let 

^x,y)^^\x-yf + /3^\xf (4.1) 
where i? and /? are two positive constants (to be chosen later), and define 
(To = sup \uo{x) - uo{y) - <^{x,y) - ^\ 

te[0,T],x,yeR'^ J 

a= sup |u(i,x) - u'(t,y) - $(a;,y) - — ^1 - (Jo, 

te[0,T],x,yeR" ^ It) 

where 7 G (0, 1). Next, we introduce the quantity 

*(t, X, y) = u{t, x) - u'{t, y) - $(x, y) - - 7^—^, (4-2) 

where 5 € (0, 1). Recall that u{t,x) and u'^{t,y) are bounded and continuous functions. Therefore, 
thanks to the penalization term there exists {to,xo,yQ) G [0,T) x M" x M" such that 

'i'{to,xo,yo) = sup 'if{t,x,y). 

t,x,y 

We are interested in finding an upper bound on cr + ctq by deriving a positive upper bound on a. 
Therefore, without loss of generality, we may assume that a > 0. This implies to > 0, since 

*(io,a;o,yo) > sup \u{t,x) - u^{t,y) ~ ^{x,y) - - 5a = ao + {1 - 6)(T > ao, 

t(E[0,T],a;,j/eR" ^ 1 — t J 

while on the other hand to = would imply that "i'{to,xo,yo) < ctq 

Therefore we can apply the maximum principle for semicontinuous functions, adapted to IPDEs[5], 
and conclude: 

For each k e (0, 1), there exists a, 6 e M satisfying 

(5(7 7 
a — = — 



T {T-tof 
such that 

a + H{to, xo, u{to, xo),'^x^{xo, yo)) < (4.3) 
b + H{to,yo,u'{to,yo),^y^{xo,yo)) +2:^($(a;o, •)) + 2:''"(u^(io, •)) > 0. (4.4) 

Note that lim„^o2:«($(a;o, •)) = 0, and u" is C^^". Therefore 

limT'''{u'{to,-)) =ei-A)-2u'{to,xo). 

Hence we pass to the limit k — > in (14. 4p . and subtract it from (|4.3p to obtain 

Sa 

T ^ (r - toY ^ ^^^o,xo,u{to,xo),V x^{xo,yo)) - H{to,yo,u''(to,yo),^y^{xo,yo)) 

-e{-A)^u'ito,xo) <0. (4.5) 

Furthermore, u{to, xo) — u'^{to, yo) ^ + T^t^j — y > as (7 > 0. Therefore, by the monotonicity 
and the Lipschitz continuity of H, we must have 



H{to,xo,u{to,xo),D^^{xo,yo)) ~ H{to,yo,u''{to,yo), Dy^{xo,yo)) > -K\xo - yo\ - Kf3\xo\. 

(4.6) 
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Furthermore, Vyu'(io,yo) = "i^lyo — xq) i.e ??|xo — yo\ < ||Vu'^||oo < C Also, it is easy to see that 
/3|xoP < C i.e l3\xo\ < Cy^- Therefore, upon denoting qo = 'd{yo — xg) and combining (|4.5p and 
we have 



^ - Xv/^ + eC{n, 1) /" [M^(to, 2/0 + z) - u\t^, yo)] < 0. (4.7) 



We now denote = e{xQ — yo) and make the substitution r = - and obtain 



-e 



[ [u%to,yo + z)- u%tQ,yo)] ^ dz 

[-""(io, 2/0 + z)- u'(to, yo) - qo-z] i^ppT 
[u'{to,yo + er) ~ u'{to,yo) - eqo.r] 



2 <e 



'S(0,1) 

= / [u^ito^yo + er) - u^^to^yo) ~ e-qo-r] . rfr 
Jb(o,i) 1^1 

where we have used the fact that u^{to,yo + z) — Ue{tQ,yo) — qo-z > — Moreover, by direct 
computation 

e\ [ [u'{to,yo + z)-u'{to,yo)]Y^dz\<Ce{l + \loge\)<Ce\\oge\ if e < e-\ 

JB{Q,e)C \Z\ 

Therefore, (j4.7l) now imphes 



Moreover, by Lipschitz continuity of uq{x) we have 

^o< sup = — -. (4.9) 

Therefore 

'-^^ < - + f + i^' + .1 log.l + C^e^. (4.10) 

2 

Now choose /? = ^/a ^2 ; and maximize the right hand side with respect to d to obtain, 

d{a + gp) 
T 

We now let (5 1 Ij and conclude 



< Ce(l + |loge|) (<Ce|loge|) if e < e' 



2 

u{t,x) - u^{t,x) ~ -^|x|2 - ^ < Ce|loge|. 
Finally, we let 7 J, in the above inequality to conclude 

u{t,x)-u\t,x) <Ce\\oge\. (4.11) 

The inequality (j4.1ip provides only half of requirement, the other half could also be concluded in 
a similar manner. 

□ 
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